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<D ' Abstract 

& 

The Vicsek model for the self-propelling particles in 2D is investigated with 

respect to addition of the stochastic perturbation of dynamic equations. We 
show that this model represents in essential the same type of bifurcations un- 
der different type of noise as the celebrated Kuramoto model of spontaneous 
synchronization. These models demonstrate similar behavior at least within 
the mean-field approach. To prove this we consider two types of noise for 
the Vicsek model which commonly considered in the literature: the intrinsic 
and the extrinsic ones (according to the terminology of J. A. Pimentel et al. 
[Phys. Rev. E 77, 061138 (2008)]). The qualitative correspondence with 
the bifurcation of stationary states in Kuramoto model is stated. New type 
£> ■ of stochastic perturbation - the "mixed" noise is proposed. It is constructed 

as the weighted superposition of the intrinsic and the extrinsic noises. The 
corresponding phase diagram "noise amplitude vs. interaction strength" is 



obtained. The possibility of the tricritical behavior for the Vicsek model is 

in 

o 
o 



predicted. 



1. Introduction 

k> ', The Vicsek model (VM) of the dynamics of self-propelling particles pro- 

posed in |IJ gave intense impact for the research of the order- disorder transi- 
tions in the systems of such kind. In this model the motion of self-propelling 
particles is driven by very simple rules. The dynamic rule proposed by Vic- 
sek is based on the constraint of alinement of particle velocity along the 
direction of the mean velocity of its neighborhood. As shown by numerous 
simulations [0, [| such dynamic rule favors the emergence of ordered motion 
at high enough densities. 
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Addition of the stochastic perturbation in the direction of motion leads 
to the transition from the ordered state of synchronized motion to the dis- 
ordered state when the intensity of the perturbation grows. As it has been 
demonstrated in j|, ||, the character of such transition depends essentially 
on the type of the stochastic perturbation. The intrinsic noise which is noth- 
ing but adding the stochastic increment to the direction of particle's motion 
originally used in [3J, leads to the continuous phase transition. The so called 
extrinsic noise introduced in [Q where some random vector is added to the 
particle velocity, makes the transition discontinuous. It should be noted that 
the same situation is well known for another model of synchronization - the 
Kuramoto model (KM) |J (see also [[F| and references therein). Yet in the 
KM the type of the transition depends on the distribution function g(u>) of 
the own frequencies of the oscillators [|7| . 

In this paper we investigate the relation between the Viscek and the 
Kuramoto models. Also we extend the results of || by introducing the novel 
type of noise perturbation which is the weighted sum of the intrinsic and the 
extrinsic ones. We call it "mixed" noise. It is shown that for such noise the 
bifurcation of stable solution becomes more complex. The tricritical behavior 



similar to that in Kuramoto model Wn is demonstrated. 



2. Relation between Vicsek and Kuramoto models 

The dynamic equation in the continuous time limit for the Vicsek model 
was considered in [|, [|. It has a simple form: 

Jt Vi = W - X V * ' (1) 

where oj v . is the "angular velocity" of z-th particle. This angular velocity 
depends on the velocities of neighboring particles. The self-propelling force 
and the frictional force are assumed to balance each other. As it is shown 
in || the angular velocity u Vt of the dynamic model, which corresponds to 
the discrete Vicsek automaton in the continuous time limit, consists of two 
parts: 

w v , = w Ui + w vu , (2) 

where u) Ui is the angular velocity corresponding to the average velocity of 
the "nearest " (in the sense of some model for the distance) neighbors of i-th 



particle: 

* (ij) 

which serves as the local order parameter. Here iVj is the number of the 
nearest neighbors of z-th particle. Another term: 

u; vu = AviX t Ui (4) 

denotes the relative angular velocity of the particle with respect to the direc- 
tion r Ui = Aof the average velocity of the neighbors @. The coefficient A 
in Eq. (§|) determines the relaxation rate (inverse of characteristic relaxation 
time) to the direction of the local order parameter ||. Such rate depends 
both on the number of the nearest neighbors and the value of the average 
velocity. Therefore it is natural that A is proportional to the local flux: 

A = \ Ui , (5) 

where A is proportional to the number of the neighbors. As it is shown in 
such dependence can be obtained in the approximation when the velocity 
and the density fields are decoupled. Throughout this paper we consider 2D 
case and use the following definition of the local order parameter: 

^ = ^E e ^'- ( 6 ) 

1 (iJ) 

In terms of this order parameter the equation (fj) takes the form || : 

e i = ip i + Asm(t/} i -e i ) , (7) 

where A is given by @. 

The dynamic equation similar to (0) appears in the theory of stochastic 
synchronization for the KM 0. The KM is determined by the following 
dynamic equations for the oscillators phases: 

9i = Ui + K ^ sin (9j - 9i) ^ Ui + K Nm sin(^-^) 5 (8) 

{hi) 

where K - is the interaction strength, ipi is the average local phase of the 
nearest oscillators |f7|. 



From the mathematical point of view dynamic equations of motion Eq. (^) 
and Eq. (^) are similar. The main difference between (|7|) and (|8]) is the form 
of the first term. For the Kuramoto model it defines the own frequency of 
the i-th oscillator. Usually the distribution of own frequencies is modeled by 
some function g{uS). For the Vicsek model the first term determines the col- 
lective contribution which is assumed to vary slowly as it is a collective mode. 
Clearly, it vanishes in the case of coherent motion, yet this term favors the 
transition to the ordered state. From this it follows that the synchronization 
in the VM is more pronounced than that in the KM with finite radius of 
interaction. Note that this term can be connected with the convective term 
v ■ Vv in the continuum approach. According to the arguments in || the 
convective term effectively generates the long-ranged correlations. In its turn 
this leads to the emergence of the ordered motion. 

Based on the mathematical similarity of the dynamic equations for the 
VM and the KM mentioned above, it is quite natural to expect that these 
models demonstrate qualitatively the same behavior at least in the mean 
field approximation. Indeed, the VM with intrinsic noise shows the super- 
critical bifurcation from disordered to ordered state |fL0 |. The same behavior 



is observed for the KM with unimodal frequency distribution of the form 
g{ui) = S(u)) [|TT|1 . Yet in the KM the type of the transition depends on the 
sign of g"(0) j7|. In the framework of the VM the distribution of the own fre- 
quencies is equivalent to the introduction of the stochastic perturbation into 
the equation of motion. Therefore for the VM there must be the dependence 
of the character of the order- disorder transition on the type of the stochastic 
perturbation. This result was obtained in numerical simulations in || and 
confirmed theoretically |5|, [L2j . 



First let us consider the case of intrinsic noise |IJ which can be modeled 
by inclusion of the Langevin source into the equation of motion with g(u) = 

d6i = -Kr sin(^) dt + \fi5dwi{t) . (9) 

Here Wi(t) is the Wiener process [O] and represents the random increment of 
the angle and D is the diffusion coefficient. In such an approach it represents 



the intensity of the intrinsic noise in the particle decision algorithm |]T2| . In 
the mean-field approximation one easily comes to the self-consistent equation 
for the stationary value of the order parameter 0: 



771^ = F (Kr) , (10) 
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Figure 1: The solution of Eq. (fio|). K = ^. Supercritical behavior of order parameter is 
shown. It is similar to the one in Vicsek model with intrinsic noise. 

where Jo and I\ are modified Bessel functions of the first kind [|14[]. This 
equation can be easily solved numerically (see Fig. [1]). The nontrivial solution 
appears continuously from the trivial one. 

As the second case we consider the extrinsic noise perturbation for the VM 
introduced in ||. This perturbation corresponds to the stochastic deviation 
of the direction of motion of i-th particle due to addition of the random 
vector ^ with the amplitude |£J = £iVj: 



6i(n + 1) = angle(wi(n) + £ £ ) 



(11) 



where £ is the control parameter. Let us define the quantity K = l/£ then the 
the self-consistent mean-field equation for the order parameter is as follows 



1 

2^ 



K r + cos a ) da 



1 + 2Kr cos a + ( Kr 



F{Kr) . 



(12) 



It is quite remarkable that ( J12) coincides exactly with corresponding equation 
for the order parameter of the Kuramoto-like model augmented with the 



phase pinning 15 



In this case the subcritical bifurcation of the solution occurs (see Fig. ^]). 
There is a discontinuous jump between ordered and disordered motion be- 
cause of the existence of the metastable state. Therefore the extrinsic noise 



F{Kr) 




Figure 2: Graphical solution of Eq. ( [12[ ) 



in the VM is equivalent to the random pinning for the direction of motion. 
For the self-propelling particle this perturbation corresponds to the instant 
change of the direction due to triggering the attention to some particle which 
moves asynchronously with the nearest neighbors forming a stable flock. Due 
to the nonlinear character of the angular dependence this "extrinsic" pertur- 
bation corresponds to the "colored" Langevin force and indeed represents 
physical situation quite different from that for the intrinsic noise 
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3. "Mixed" noise 

The similarity between the bifurcations of the stationary states in the 
Vicsek and the Kuramoto models is based on the equivalence of their equa- 
tions of motion. The fact that the Vicsek model with the extrinsic noise 
demonstrates the hysteretic behavior typical for the first order phase transi- 
tion and continuum phase transition for the intrinsic noise, allows us to con- 
clude that these perturbations can be present simultaneously in the general 



case of the dynamics. Therefore the tricritical behavior might be expected if 
some switching parameter between these two perturbations is added. Note 
that tricritical behavior is known for the Kuramoto model with symmetrical 
bimodal distribution function 

g v (ui) = r)8(u> - u ) + (1 - T])6(u + u> ), (13) 

with rj = 1/2 (see [0). 

The simplest variant is the statistical mixing similar to that of (|13|) the 
extrinsic and the intrinsic noises. The corresponding self-consistent equation 
is of the form: 

r = 77 F v (r) + (1 - rj) F,(r) . (14) 

Here the right hand side of (|HD corresponds to the statistically mixed prob- 
ability density for the angle of direction: 

fmixie, r) = r)f v (e, r) + (i - v )f s (e, r) , (15) 

where f s and /„ are distribution functions for intrinsic and extrinsic noises 
respectively and < rj < 1 denotes the parameter of statistical mixing. 

As the analysis of ([TJ|) shows, the bifurcation of its solutions with varying 
parameter rj corresponding to the tricritical behavior (see Fig. ^j. Note that 
the same behavior is obtained for the KM if the number of modes included 



into frequency distribution function g(u) more than one (see [|16], pj 

Considering ( JT5| ) as the stationary solution of the Fokker-Plank equation 
and using the standard relation with the Langevin equation (see e.g. [p"3"[): 



d9i = -K r sm{8i) dt + b{6i,r)dwi{t) (16) 

one can show that: 

e 

b 2 (6,r;r]) = — — Kr f mix (x,r) sinxdx , (17) 

where 8 can be found from the condition that in accordance with (|9]) for the 
intrinsic noise (77 = 0) b(6, r) = 1 : 

/In 2 

9 = arccos — - — h cos 9 
\ Kr 




2.20 



Figure 3: Solutions for the order parameter r vs. K. Gray region shows the situation 
when two stable regimes are possible. The value of mixing parameter 77 rs 0.16. 

Then from ([T7|) one can easily get the following expression for b(9,r; 77): 

b 2 (6, r; rf) = (1 - 77) + v(B ex t(0) ~ B ext (6 )) , (18) 

where the function: 



Bext(Q) 



2r 2 / Vl + r 2 + 2rcos# , ^ 

arctan I , | — yr 1 + 2r cos# + 1 , 



VT+ 






VTT 






determines the corresponding parameter of Fokker-Plank equation in case of 
extrinsic noise (77 = 1). 

As the analysis shows, the addition of "mixed" noise leads to complex 
behavior. It is natural that for 77 — > 1 the behavior corresponds to one with 
intrinsic noise where the supercritical bifurcation takes place. Small values of 
mixing parameter 77 drives the system into the situation with intrinsic noise 
for which the subcritical regime occurs. By appropriate tuning the parameter 
77 the possibility when both behaviors supercritical and subcritical exist at one 
time can be realized. Indeed for the case of "mixed" noise there is a region of 
values K and 77 where two nontrivial stable regimes exist simultaneously (see 
Fig. [|). The complete phase diagram for the mixed noise is shown in Fig. |5|. 
Line AO stands for the discontinuous phase transition (like in extrinsic noise 
case). Line BC is for the continuous phase transition (like in case of the 
intrinsic noise). The point O is the tricritical point, when lines of phase 
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Figure 4: Graphical solution of self-consistent equation for the case of mixed noise with 
different values of parameters r\ and K 
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Figure 5: Phase diagram for "mixed" noise. 



transitions of different orders meet (K = 2.0, 77 ~ 0.28). Thus the statistical 
mixing of the extrinsic and the intrinsic noises leads to the emergence of the 
tricritical point. This prediction can be tested in computer simulation of the 
corresponding dynamics. 

Note that the function g(u) plays a role of initial distribution for the 
frequencies and represents the background noise in the KM. In the Vicsek 
model there is no direct analog of such function. Yet as has been shown 
above, there is a correspondence between the bifurcation behavior of the KM 
determined by g(cu) (unimodal, bimodal, etc.) and the VM with specific 
kind of noise. It is interesting to find the relation between the properties of 
g(ui) and the type of noise in the VM which leads to the similarity of the 
phase diagram. The question whether the similarity between these models 
conserves with the inclusion of fluctuations, remains open. 

4. Discussion 

The results of the work can be summarized as follows. It is shown that 
the type of ordering for the Vicsek model depends on the kind of noise per- 
turbation of the equation of motion. The transition to ordered state can be 
either of continuous or discontinuous for intrinsic and extrinsic noises cor- 
respondingly. It is demonstrated that the Vicsek model shows qualitatively 
the same variety of the phase behavior as the Kuramoto model of the phase 
synchronization of nonlinearly coupled oscillators. Both known cases of ex- 
trinsic and intrinsic noise action in Vicsek model demonstrate the results 
similar to those for Kuramoto model. In particular on the basis of the stated 
equivalence of these models "mixed" noise is investigated and the tricritical 
behavior for the Vicsek model is demonstrated. 
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